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Abstract 



We investigate the flow of two immiscible, viscous fluids in a rotating Hele- 
Shaw cell, when one of the fluids is a ferrofluid and an external magnetic field is 
applied. The interplay between centrifugal and magnetic forces in determining 
the instability of the fluid-fluid interface is analyzed. The linear stability 
analysis of the problem shows that a non-uniform, azimuthal magnetic field, 
applied tangential to the cell, tends to stabilize the interface. We verify that 
maximum growth rate selection of initial patterns is influenced by the applied 
field, which tends to decrease the number of interface ripples. We contrast 
these results with the situation in which a uniform magnetic field is applied 
normally to the plane defined by the rotating Hele-Shaw cell. 
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When a fluid is pushed by a less viscous one in a narrow space between two parallel 
plates (the so-called Hele-Shaw cell) , the well-known Saff man- Taylor instability phenomenon 
arises PJ, which takes the form of fingering 0. Traditionally, experiments and theory focus 
on two basic Hele-Shaw flow geometries: (i) rectangular |J and radial ||. In rectangular 
geometry cells the less viscous fluid is pumped against the more viscous one along the 
direction of the flow. Meanwhile, in the radial geometry case, the less viscous fluid is 
injected to invade the more viscous one, through an inlet located on the top glass plate. In 
both geometries, the viscosity- driven instability leads to the formation of beautiful fingering 
patterns. 

In recent years, the quest for new morphologies and richer dynamic behavior resulted in 
a number of modifications of the classic Saffman- Taylor setup [[|]. An interesting variation 
of the traditional viscosity-driven fingering instability is the investigation of radial Hele- 
Shaw flows in the presence of centrifugal driving. The inclusion of centrifugal forces can be 
considered by rotating the cell, with constant angular velocity, around an axis perpendicular 
to the plane of the flow. In this case, the interface instability can be driven by the density 
difference between the fluids. In the late 1980s Schwartz |5j performed the linear stability 
analysis of the rotating cell problem, in the limits of high density and viscosity contrast. 
More recently, Carrillo et al. || studied, both theoretically and experimentally, flow in a 
rotating Hele-Shaw cell in arbitrary density and viscosity contrast. They extended the linear 
analysis performed in || by considering that the inner fluid is injected in the cell through a 
hole at the center of rotation, with constant injection rate. The linear growth rate calculated 
in || shows that the interface instability can be driven by both the density difference and the 
viscosity contrast between the fluids. Their experimental results supported their theoretical 
analysis. Carrillo and co-workers also examined the radial displacement of a rotating fluid 
annulus, bound by a second fluid, in stable |7[ and unstable || regimes. In another interesting 
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work, Magdaleno et al. || applied a conformal mapping technique to argue that the effects 
of rotation can be used to prevent cusp singularities in zero surface tension Hele-Shaw flows. 

Another stimulating modification of the traditional Saffman- Taylor problem in non- 
rotating Hele-Shaw cells considers the interface morphology when one of the fluids is a 



ferrofluid |]10[ , and an external magnetic field is applied perpendicular to the cell plates. 
Ferrofluids, which are colloidal suspensions of microscopic permanent magnets, respond 
paramagnetically to applied fields. As a result of the ferrofluid interaction with the external 
field, the usual viscous fingering instability is supplemented by a magnetic fluid instabil- 
ity [pro | , resulting in a variety of new interfacial behaviors [|IT| -|T6|j. 

The richness of new behaviors introduced by both rotation and magnetic field into the 
traditional Saffman- Taylor problem, motivated us to analyze the situation in which these 
two effects are simultaneously present. In this work we perform the linear stability analysis 
for flow in a rotating Hele-Shaw cell, assuming that one of the fluids is a ferrofluid and that 
a magnetic field is applied. First, we consider the new situation in which a non-uniform, 
azimuthal, in-plane field is applied. The competition between rotation and magnetic field is 
analyzed. We show the azimuthal magnetic field provides a new mechanism for stabilizing 
the interface. It is the first time such a field induced, stabilizing mechanism is proposed, 
and proved to be effective for flow in rotating Hele-Shaw cells. Through the analysis of 
the maximum growth rate, we verify that the azimuthal magnetic field acts to decrease the 
number of interface ripples. Finally, we contrast these results with the destabilizing situation 
in which a uniform magnetic field is applied normally to the plane defined by the rotating 
Hele-Shaw cell. 

Consider a Hele-Shaw cell of thickness b containing two immiscible, incompressible, vis- 
cous fluids (see figure 1). Denote the densities and viscosities of the inner and outer fluids, 
respectively as pi, T]i and p 2 , r] 2 . The flows in fluids 1 and 2 are assumed to be irrotational. 
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Between the two fluids there exists a surface tension a . We assume that the inner fluid is the 
ferrofluid (magnetization M), while the outer fluid is non-magnetic. During the flow, the 
fluid-fluid interface has a perturbed shape described as 1Z = R + ((9, t), where 9 represents 
the polar angle, and R is the radius of the initially unperturbed interface. 

In order to include centrifugal forces, we allow the cell to rotate, with constant angular 
velocity Q, about an axis perpendicular to the plane of the flow (figure 1). To include 
magnetic forces, we consider the action of an external magnetic field H, produced by a long, 
straight wire carrying a current J, directed along the axis of rotation. By Ampere's law, it 
can be shown that the steady current I produces an azimuthal magnetic field external to 
the wire H = (I /2irr) 9, where r is the distance from the wire, and 9 is the unit vector 
pointing in the direction of increase of 9. 



Following the standard approximations used by Rosensweig Jl0| and others |]Tl"|-|T3|j we 
assume that the ferrofluid magnetization M is collinear with the external field H and that 
the influence of the demagnetizing field is neglected. It is also assumed that the ferrofluid is 
electrically nonconducting and that the displacement current is negligible. For the quasi two- 
dimensional geometry of a Hele-Shaw cell, the three dimensional flow may be replaced with 
an equivalent two-dimensional flow v(x, y) by averaging over the z direction perpendicular 
to the plane of the Hele-Shaw cell. Imposing no-slip boundary conditions and a parabolic 



velocity profile one derives Darcy's law for ferrofluids in a Hele-Shaw cell []T3 , |17 |, which must 
be augmented by including centrifugal forces, 

b 2 { - l r+ b / 2 



r)v = |Vp - ^ jT^ no(M ■ V)Hdz - P n 2 r f j , (1) 

where p is the hydrodynamic pressure, /io is the free-space permeability, and f denotes a 
unit vector pointing radially outward. Equation ([[]) describes non-magnetic fluids by simply 
dropping the terms involving magnetization. 

Since the velocity field v is irrotational, it is convenient to rewrite equation ([[]) in terms 
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of velocity potentials. We write v = — V0, where defines the velocity potential. Similarly, 
we rewrite the magnetic body force in equation (fl|) as Ho{M ■ V)H = [IqMVH = V^, where 
we have introduced the scalar potential 

y = H jM(H)dH=^^, (2) 

with M = M(H) = xH, % being a constant magnetic susceptibility. 

With the definitions of and \l/ we notice that both sides of equation (H) are recognized 
as gradients of scalar fields. After integrating both sides of equation (|T|), we evaluate it for 
each of the fluids on the interface. Then, we subtract the resulting equations from each 
other, and divide by the sum of the two fluids' viscosities to get the equation of motion 

To obtain (|3|) we have used the pressure boundary condition pi — p 2 = ctk at the interface, 
where k = {[r 2 + 2{dr/d9f - r(d 2 r/de 2 )}/[r 2 + {dr/deff 2 } denotes the interfacial curva- 
ture in the plane of the Hele-Shaw cell. The dimensionless parameter A = 7 ?i)/( 7 ?2 + ^i) 
is the viscosity contrast. 

For the purpose of the following linear analysis, we perturb the interface with a single 
Fourier mode 

C(9,t) = Cn(t)exp(in9), n = 0,l,2,... (4) 

The velocity potential for fluid j (j = 1,2 indexes the inner and outer fluids, respectively), 
<f)j, obeys Laplace's equation V 2 0j = and can be written as 

<f>j = 4>°j+ <Pjn i—j exp(m^), (5) 

where 0° are independent of r and 6. 

We need additional relations expressing the velocity potentials <pj in terms of the per- 
turbation amplitudes (, in order to conclude our derivation and close equation @. To find 
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these, we considered the kinematic boundary condition |T(J, which refers to the continuity 



of the normal velocity across the interface. Inserting expression (|J) for ((0,t) and (|5|) for cf)j 
into the kinematic boundary condition, we solved for <pj n consistently to first order in ( to 
find 

cj> jn = (-l)'-Cn, (6) 

n 

where the overdot denotes total time derivative. 

Substitute expression (^j) for <pj n into equation of motion ([|), and again keep only linear 
terms in the perturbation amplitude. This procedure eliminates the velocity potentials 
from equation (0), and we obtain the differential equation for the perturbation amplitudes 
( n = X(n)( n , implying that the relaxation or growth of the mode n is proportional to the 
factor exp[A(n)t], where 

A(n) = 12 i N ° - Nb ~ ( " 2 - i} ] < 7 > 

is the linear growth rate. We define the dimensionless parameters Nq = [-R 3 (pi — p2)Q 2 ]/o', 
and A^b = A i oX^ 2 /(4vr 2 cri?) as the rotational and magnetic bond numbers, respectively. Nq 
{Nb) measures the relative strength of centrifugal (magnetic) and capillary effects. 

Inspecting equation (0) for the linear growth rate A(n) we observe the interplay of rota- 
tion, magnetic field and surface tension in determining the interface instability. If A(n) > 
the disturbance grows, indicating instability. As usual, the contribution coming from the 
surface tension term has a stabilizing nature (a stabilizes modes of large n). The factor 
(n 2 — 1) in equation (0) arises directly from the first order terms in ( present in the curva- 
ture, while the overall factor of n can be traced to the fact that in the generalized Darcy's 
law (0) the velocity is proportional to gradients of an effective pressure. With these con- 
siderations in mind, let us focus on the relation between rotation and magnetic field. As a 
result of centrifugal forcing Nq may be either positive or negative, depending on the rela- 
tive values of the fluid's densities. If the inner fluid is more dense (pi > p 2 ), Nq > and 
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rotation plays a destabilizing role. The opposite effect arises when p\ < pi. On the other 
hand, the azimuthal magnetic field contribution Nb always tends to stabilize the interface. 
This indicates that the rotation-driven instability could be delayed or even prevented if a 
sufficiently strong non-uniform, azimuthal magnetic field is applied in the plane of the flow. 

A physical explanation for the stabilizing role of the magnetic field can be given based 
on its symmetry properties and non-uniform character. Notice that such a field possesses a 
radial gradient. The magnetic field influence is manifested as the existence of a body force 
due to field nonuniformity. The field produces a force directed radially inward, that tends to 
move the ferrofluid toward the current- carrying wire (regions of higher magnetic field). This 
force opposes the centrifugal force and favors interface stabilization. This effect is similar 



to the gradient-field stabilization mechanism discussed by Rosensweig [10] for inviscid three 
dimensional fluid flow problems, and by Zahn and Rosensweig [|18[ for viscous, unconfined 
ferrofluids. 

In order to illustrate the role of the magnetic effects in the linear stages of the interface 
evolution, we plot in figure 2 the dimensionless growth rate X(n) = [12(r/i + r]2)R 3 /b 2 <j] X(n) 
as a function of the mode number n, taking Nq=200 and Nb=0, 100, 200. As expected, 
for zero applied field (Nb = 0) the interface is unstable. If we increase the magnitude of 
the applied field (Nb = 100), even though we get a narrower band of unstable modes, the 
interface remains unstable. If we keep increasing the magnetic field intensity, A(n) becomes 
negative for all n > 1, and the interface tends to stabilize. In figure 2 we see that for 
Nb = 200 the interface is stable, due to the action of the azimuthal magnetic field. 

A relevant physical quantity can be extracted from the linear growth rate: the fastest 
growing mode n*, given by the closest integer to the maximum of equation (^) with respect 
to n (defined by setting d\{n)/dn = 0), 



n max = Jl[l + (N n -N B )}. (8) 



As discussed by Carrillo et al. || for the non- magnetic, rotating case, n* is strongly correlated 
to the number of ripples present in the early stages of pattern formation. With the help of 
their experiments, the authors in || compared the number of ripples with n* and found a 
remarkable agreement. Taking this fact into account, and inspecting equation (|8]) we verify 
that, for positive Nq, an increasingly larger Ng does not only tend to decrease the finger 
growth rate, but also tends to decrease the number of interface ripples. The azimuthal 
magnetic field (or correspondingly, Nb) can be seen as a control parameter to discipline the 
number of interface undulations. 

We conclude by contrasting the results obtained above (non-uniform, tangential field) 
with those which arise when a uniform external magnetic field is applied perpendicular to a 
rotating Hele-Shaw cell containing ferrofluid. By performing the linear stability analysis of 
the system, and using the closed form expressions for the magnetic term, recently derived 



in reference |[14|| , we obtain the following linear growth rate 

b 2 an 
12(77! +772)^3 



fy^ 1 (TTI r 1 

A± H = 7777 I , m Nn + V n {p)N^ - (n 2 - 1) , (9) 



where 



Vnip) = y 



ip [n H — — ip 



1) \2 



+ 



Qn-1/2 \ P t J - Ql/2 \ -t 



(10) 



p I \ p 

Ng = fi M 2 b/27ra is the magnetic bond number for the perpendicular field configuration, 
and p = 2R/b is the aspect ratio. The aspect ratio p should not be confused with the pressure. 
Q n represents the Legendre function of the second kind, while the Euler's psi function ip is 
the logarithmic derivative of the Gamma function JOJ. Notice that the function T> n (p) > 
for n > 0. In contrast to the non-uniform, azimuthal applied field discussed earlier, a 
uniform, perpendicular magnetic field tends to destabilize the interface. If the interface is 
already unstable with respect to rotations (N n > 0), the introduction of the magnetic field 
increases the interface instability even further. On the other hand, if the outer fluid is more 
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dense (Nq < 0), the interface can still be deformed by the action of a perpendicular magnetic 
field. 

By comparing the perpendicular field growth rate expression ([)]) with its azimuthal 
field counterpart (|7|), we notice that the presence of the function T> n {p) in (|]) increases 
the complexity of the problem. Simple expressions for the azimuthal case may become 
not so simple in the perpendicular field situation. Consider, for example, the calculation 
of the fastest growing mode n*. In contrast to the simple expression for n max obtained 
in the azimuthal case (see equation (|8])), the equivalent expression for the perpendicular 
configuration cannot be written in a simple form. It is now given by the solution of the 
transcendental equation 

1 



n 2 



l+N a + ^-(n V n (pj)Ni 



(11) 



dn 

Numerical evaluation of equation ( |TTD shows that n max (and consequently, n*) increases with 
the magnitude of the perpendicular applied field. In this case, larger values of Ng tend to 
increase the number of interface ripples (see figure 3). 

In summary, we have shown that the inclusion of magnetic effects into the rotating cell 
problem provides new mechanisms for stabilizing/destabilizing the interface. Based on the 
relative simplicity of the experimental Hele-Shaw setup, it would be of considerable interest 
to perform experiments in Hele-Shaw cells simultaneously including centrifugal, magnetic 
and injection driving. The combination of these effects is likely to lead to new and exciting 
interfacial patterns in the highly nonlinear regime. 

This work was supported by CNPq and FINEP. 
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FIGURE CAPTIONS 



FIG. V. Schematic configuration of the rotational Hele-Shaw flow with ferrofmid. 

FIG. 2: Variation of the dimensionless growth rate A(n) = [12 (771 + r]2)R 3 /b 2 a] A(n) as a 
function of n for N n = 200 and (a) N B = 0, (b) N B = 100, and (c) N B = 200. The peak 
location and width of the band of unstable modes decrease with increasing 7V#. 

FIG. 3: Plot of n max as a function of iV^- for increasing values of Nn, and p = 20. These 
curves are obtained by numerically solving equation (|TT|). 
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